function [dPsiM] = dPsiPNbydM(f, Params)
% 
% DPSIPNBYM - compute the derivative of the 3.5PN phase with respect
% to the total mass M. 
% 
% usage: [dPsiM] = dPsiPNbydM(f, Params)
%
% P. Ajith, 05.04.2008
%
% $Id: dPsiPNbydM.m 59 2010-01-20 21:41:01Z anand.sengupta $

    MSOLAR_TIME = 4.92579497077314e-06;

    m = Params.totalMass.*MSOLAR_TIME;
    eta = Params.eta;

dPsiM = (-37.8569749799331.*f)./power(f.*m,0.6666666666666666) - ...
  55./(384.*f.*power(m,2).*pi) - (27145.*f)./ ...
   (64512.*power(f.*m,1.3333333333333333).*power(pi,0.3333333333333333)) -  ...
  (65.*pi)./(384.*m) + (378515.*f.*power(pi,1.6666666666666667))./ ...
   (96768.*power(f.*m,0.3333333333333333)) -  ...
  (3.183802571673538.*f)./(power(f.*m,0.6666666666666666).*eta) -  ...
  (5.*f)./(128.*power(f.*m,2.6666666666666665).*power(pi,1.6666666666666667).* ...
     eta) - 3715./(32256.*f.*power(m,2).*pi.*eta) -  ...
  (15293365.*f)./(6.5028096e7.*power(f.*m,1.3333333333333333).* ...
     power(pi,0.3333333333333333).*eta) +  ...
  (f.*power(pi,0.3333333333333333))./ ...
   (4.*power(f.*m,1.6666666666666667).*eta) -  ...
  (107.*power(f.*m,0.3333333333333333).*power(pi,0.3333333333333333))./ ...
   (42.*m.*eta) + (38645.*pi)./(32256.*m.*eta) +  ...
  (77096675.*f.*power(pi,1.6666666666666667))./ ...
   (1.6257024e7.*power(f.*m,0.3333333333333333).*eta) -  ...
  (3085.*f.*eta)./ ...
   (9216.*power(f.*m,1.3333333333333333).*power(pi,0.3333333333333333)) +  ...
  (76055.*f.*power(pi,0.3333333333333333).*eta)./ ...
   (221184.*power(f.*m,0.6666666666666666)) -  ...
  (74045.*f.*power(pi,1.6666666666666667).*eta)./ ...
   (48384.*power(f.*m,0.3333333333333333)) -  ...
  (127825.*f.*power(pi,0.3333333333333333).*power(eta,2))./ ...
   (165888.*power(f.*m,0.6666666666666666)) -  ...
  (107.*f.*power(pi,0.3333333333333333).* ...
     log(power(f.*m,0.3333333333333333).*power(pi,0.3333333333333333)))./ ...
   (42.*power(f.*m,0.6666666666666666).*eta);

